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Abstract. We show that, in a model of modified gravity based on the 
spectral action functional, there is a nontrivial coupling between cosmic 
topology and inflation, in the sense that the shape of the possible slow- 
roll inflation potentials obtained in the model from the nonperturbative 
form of the spectral action are sensitive not only to the geometry (flat 
or positively curved) of the universe, but also to the different possible 
non-simply connected topologies. We show this by explicitly comput- 
ing the nonperturbative spectral action for some candidate flat cosmic 
topologies given by Bieberbach manifolds and showing that the resulting 
inflation potential differs from that of the flat torus by a multiplicative 
factor, similarly to what happens in the case of the spectral action of 
the spherical forms in relation to the case of the 3-sphere. We then show 
that, while the slow-roll parameters differ between the spherical and flat 
manifolds but do not distinguish different topologies within each class, 
the power spectra detect the different scalings of the slow-roll potential 
and therefore distinguish between the various topologies, both in the 
spherical and in the flat case. 



1. Introduction 

Noncommutative cosmology is a new and rapidly developing area of re- 
search, which aims at building cosmological models based on a "modified 
gravity" action functional which arises naturally in the context of noncom- 
mutative geometry, the spectral action functional of |4] . As we discuss more 
in detail in ^ below, this functional recovers the usual Einstein-Hilbert 
action, with additional terms, such as a conformal gravity, Weyl curvature 
term. It also has the advantage of allowing for interesting couplings of 
gravity to matter, when extended from manifolds to "almost commutative 
geometries" as in 17J and later models |6], [3]. Thus, this approach makes 
it possible to recover from the same spectral action functional, in addition 
to the gravitational terms, the full Lagrangian of various particle physics 
models, ranging from the Minimal Standard Model of |7j, to the extension 
with right handed neutrinos and Majorana mass terms of [6J, and to su- 
persymmetric QCD as in p]. The study of cosmological models derived 
from the spectral action gave rise to early universe models as in jT5] and 
[TT] . which present various possible inflation scenarios, as well as effects on 
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primordial black holes evaporation and gravitational wave propagation. Ef- 
fects on gravitational waves, as well as inflation scenarios coming from the 
spectral action functional, were also recently studied in [T7j, [IH], 

Our previous work [TB] showed that, when one considers the nonpertur- 
bative form of the spectral action, as in [S], one obtains a slow-roll potential 
for inflation. We compared some of the more likely candidates for cosmic 
topologies (the quaternionic and dodecahedral cosmology, and the flat tori) 
and we showed that, in the spherical cases (quaternionic and dodecahedral), 
the nonperturbative spectral action is just a multiple of the spectral action 
of the sphere S^, and consequently the inflation potential only differs from 
the one of the sphere case by a constant scaling factor, which cancels out in 
the computation of the slow-roll parameters, which are therefore the same 
as in the case of a simply connected topology and do not distinguish the 
different cosmic topologies with the same spherical geometry. 

This result for spherical space forms was further confirmed and extended 
in [25] , where the nonperturbative spectral action is computed explicitly for 
all the spherical space forms and it is shown to be always a multiple of the 
spectral action of S^, with a proportionality factor that depends explicitly on 
the 3-manifold. Thus, different candidate cosmic topologies with the same 
positively curved geometry yield the same values of the slow-roll parameters 
and of the power-law indices and tensor-to-scalar ratio, which are computed 
from these parameters. 

In jl6j . however, we showed that the inflation potential obtained from 
the nonperturbative spectral action is different in the case of the flat tori, 
and not just by a scalar dilation factor. Thus, we know already that the 
possible inflation scenarios in noncommutative cosmology depend on the 
underlying geometry (flat or positively curved) of the universe, and the 
slow-roll parameters are different for these two classes. 

The slow-roll parameters alone only distinguish, in our model, between 
the flat and spherical geometries but not between different topologies within 
each class. However, in the present paper we show that, when one looks at 
the amplitudes for the power spectra for density perturbations and gravi- 
tational waves (scalar and tensor perturbations), these detect the different 
scaling factors in the slow-roll potentials we obtain for the different spherical 
and flat topologies, hence we obtain genuinely different inflation scenarios 
for different cosmic topologies. 

We achieve this result by relying on the computations of the nonpertur- 
bative spectral action, which in the spherical cases are obtained in |16| and 
[25], and by deriving in this paper the analogous explicit computation of the 
nonperturbative spectral action for the flat Bieberbach manifolds. A similar 
computation of the spectral action for Bieberbach manifolds was simultane- 
ously independently obtained by Piotr Olczykowski and Andrzej Sitarz in 

m- 

Thus, the main conclusion of this paper is that a modified gravity model 
based on the spectral action functional predicts a coupling between cosmic 
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topology and inflation potential, with different scalings in the power spectra 
that distinguish between different topologies, and slow-roll parameters that 
distinguish between the spherical and flat cases. 

The paper is organized as follows. We first describe in the broader 
context in which the problem we consider here falls, namely the cosmological 
results relating inflation, the geometry of the universe, and the background 
radiation, and the problem of cosmic topology. We then review briefly in ^ 
the use of the spectral action as a modified gravity functional and the impor- 
tant distinction between its asymptotic expansion at large energies and the 
nonperturbative form given in terms of Dirac spectra. In ^ we present the 
main mathematical result of this paper, which gives an explicit calculation 
of the nonperturbative spectral action for certain Bieberbach manifolds, us- 
ing the Dirac spectra of [21] and a Poisson summation technique similar to 
that introduced in [5], and used in [T6] and [25]. 

Finally, in ^we compare the resulting slow-roll inflation potentials, power 
spectra for density perturbations and slow-roll parameters, for all the dif- 
ferent possible cosmic topologies. 

1.1. Inflation, geometry, and topology. It is well known that the mech- 
anism of cosmic inflation, first proposed by Alan Guth and Andrei Linde, 
naturally leads to a flat or almost flat geometry of the universe (see for 
instance §1.7 of |14|). It was then shown in [10] that the geometry of the 
universe can be read in the cosmic microwave background radiation (CMB), 
by showing that the anisotropies of the CMB depend primarily upon the 
geometry of the universe (flat, positively or negatively curved) and that this 
information can be detected through the fact that the location of the first 
Doppler peak changes for different values of the curvature and is largely 
unaffected by other parameters. This theoretical result made it possible to 
devise an observational test that could confirm the inflationary theory and 
its prediction for a flat or nearly flat geometry. The experimental confir- 
mation of the nearly flat geometry of the universe came in [2j through the 
Boomerang experiment. Thus, the geometry of the universe leaves a mea- 
surable trace in the CMB, and measurements confirmed the flat geometry 
predicted by inflationary models. 

The cosmic topology problem instead concentrates not on the question 
about the curvature and the geometry of the universe, but on the possible 
existence, for a given geometry, of a non-simply connected topology, that 
is, of whether the spatial sections of spacetime can be compact 3-manifolds 
which are either quotients of the 3-sphere (spherical space forms) in the 
positively curved case, quotients of 3-dimensional Euclidean space (flat tori 
or Bieberbach manifolds) in the flat case, or quotients of the 3-dimensional 
hyperbolic space (hyperbolic 3-manifolds) in the negatively curved space. A 
general introduction to the problem of cosmic topology is given in [12j . 

Since the cosmological observations prefer a flat or nearly flat positively 
curved geometry to a nearly flat negatively curved geometry (see |2], [27]). 
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most of the work in trying to identify the most likely candidates for a non- 
trivial cosmic topology concentrate on the flat spaces and the spherical space 
forms. Various methods have been devised to try to detect signatures of cos- 
mic topology in the CMB, in particular through a detailed analysis of simu- 
lated CMB skies for various candidate cosmic topologies (see [22J for the flat 
cases). It is believed that perhaps some puzzling features of the CMB such 
as the very low quadrupole, the very planar octupole, and the quadrupole- 
octupole alignment may find an explanation in the possible presence of a 
non-simply connected topology, but no conclusive results to that effect have 
yet been obtained. 

The recent results of [16] show that a modified gravity model based on the 
spectral action functional imposes constraints on the form of the possible 
inflation slow-roll potentials, which depend on the geometry and topology of 
the universe, as shown in [16j. While the resulting slow-roll parameters and 
spectral index and tensor-to-scalar ratio distinguish the even very slightly 
positively curved case from the flat case, these parameters alone do not 
distinguish between the different spherical topologies, as shown in [25]. As 
we show in this paper, the situation is similar for the flat manifolds: these 
same parameters alone do not distinguish between the various Bieberbach 
manifolds (quotients of the flat torus), but they do distinguish these from 
the spherical quotients. However, if one considers, in addition to the slow- 
roll parameters, also the power spectra for the density fluctuations, one can 
see that, in our model based on the spectral action as a modified gravity 
functional, the resulting slow-roll potentials give different power spectra that 
distinguish between all the different topologies. 

1.2. Slow-roll potential and power spectra of fluctuations. We first 
need to recall here some well known facts about slow-roll inflation potentials, 
slow-roll parameters, and the power spectra for density perturbations and 
gravitational waves. We refer the reader to [23] and to |13j . |24j . as well as 
to the survey of inflationary cosmology fT]. 

Consider an expanding universe, which is topologically a cylinder y x M, 
for a 3-manifold Y, with a Lorentzian metric of the usual Priedmann form 

(1.1) ds^ = -dt^ + a{t)'^dsl 

where ds\ = gijdx^dx^ is the Riemannian metric on the 3-manifold Y. 
In models of inflation based on a single scalar field slow-roll potential 



y(0), the dynamics of the scale factor a{t) in the Friedmann metric (1.1) is 
related to the scalar field dynamics through the acceleration equation 

(1.2) ^ = H\l-e), 

a 

where H is the Hubble parameter, which is related to the scalar field and 
the inflation potential by 

(1.3) = ^(^^4>^ + V{^)y and 4> + 3H4> + V {((>) = 0, 
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and e is the slow-roll parameter, which depends on the potential V as de- 



scribed in (1.11) below, see pTj for more details. 



It is customary to decompose perturbations of the metric ds of (1.1) into 
scalar and tensor perturbations, which correspond, respectively, to density 
fluctuations and gravitational waves. One typically neglects the remaining 
vector components of the perturbation, assuming that these are not gener- 
ated by inflation and decay with the expansion of the universe, see §9.2 of 
[1]. Thus, one writes scalar and tensor perturbations in the form 

(1.4) ds^ = -{l+2^)dt^+2a{t)dBdt+a{tf{{l-2^)gij+2/\E+hij)dx^dx^ 

with dB = diBdx^ and Ai? = didjE, and where the hij give the tensor 
part of the perturbation, satisfying d^hij = and = 0. The tensor 
perturbations hij have two polarization modes, which correspond to the two 
polarizations of the gravitational waves. 

One considers then the intrinsic curvature perturbation 

(1.5) 11 = ^-^6^, 

which measures the spatial curvature of a comoving hypersurface, that is, a 
hypersurface with constant <j). After expanding TZ in Fourier modes in the 
form 

(1.6) K = 

one obtains the power spectrum Vs{k) for the density fluctuations (scalar 
perturbations of the metric) from the two-point correlation function, 

(1.7) {UkUk') = (27r2)3 Vs{k) 5\k + k'). 

In the case of a Gaussian distribution, the power spectrum describes the 
complete statistical information on the perturbations, while the higher or- 
der correlations functions contain the information on the possible presence 
of non-Gaussianity phenomena. The power spectrum Vt{k) for the tensor 
perturbations is similarly obtained by expanding the tensor fluctiations in 
Fourier modes hk and computing the two-point correlation function 

(1.8) {hkhk') = (2^2)3 p^(^) ^3(^ ^ i^fy 

See [T], §9.3, and [24J for more details. 

In slow-roll inflation models, the power spectra Vs{k) and Vt{k) are re- 
lated to the slow-roll potential V{(j)) through the leading order expression 
(see [,23j) 

(1-9) 'Psik)^^^ and Vt{k) ^ 



M%{V'Y ' Ml,,' 

up to a constant proportionality factor, and with Mpi the Planck mass. Here 
the potential V{4') and its derivative V'{(j)) are to be evaluated at A; = aH 
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where the corresponding scale leaves the horizon during inflation. These can 
be expressed a power law as ([23]) 

-ns + ^log(k/ko) 



(1.10) 



Vtik)^ Vtiko) 
where the spectral parameters Hg, 



rif, ttsi and at depend on the slow-roll 
potential in the following way. In the slow-roll approximation, the slow-roll 
parameters are given by the expressions 

v 



lOvr 



(1.11) 



Mil y" 



V 

V'V' 



(1.12) 



nt 



as 



2C 



at ~ 



^ 647r2 y2 

Notice that we follow here a different convention with respect to the one we 
used in [16] on the form of the slow-roll parameters. The spectral parameters 
are then obtained from these as 

l-Qe + 2r] 
-2e 

16e?7 - 24e2 
4er? - 8e2 

while the tensor-to-scalar ratio is given by 

Ps 

From the point of view of our model, the following observation will be 
useful when we compare the slow-roll potentials that we obtain for different 
cosmic topologies and how they affect the power spectra. 

Lemma 1.1. Given a slow-roll potential V{(f)) and the corresponding power 
spectra Vs{k) and Vt{k) as in (1.9) and ( |1.10| ). // the potential V{(t)) is 
resettled by a constant factor V{<j)) i— )• Ay(0), then the power spectrtt Vs{k) 
and Vt{k) are also resettled by the same factor A > 0, while in the power law 
( 1.10| ) the exponents are unchttnged. 



(1.13) 



16e. 



Proof. This is an immediate consequence of (1.9), ( 1.11[ ), ( |1.12 ), and (1.10). 
In fact, from (1.9), we see that V^cp) i— )• XV{(j)) maps Vt ^ XPt, and also 
Vs ^ X'Ps, since it transforms V'^CV)"'^ ^ XV'^iy')^'^. On the other hand, 
the expressions {V /V)'^ and V" /V and V'V" /V'^ in the slow-roll parameters 



(1.11) are left unchanged hy V XV, so that the slow-roll parameters and 



all the resulting spectral parameters of (1.12) are unchanged. Thus, the 
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power law (1.10) only changes by a multiplicative factor Vs{ko) i— )• XVs{ko) 



and Vt{ko) i— )• XVt{ko), with unchanged exponents. □ 

2. NONCOMMUTATIVE COSMOLOGY 

2.1. The spectral action as a modified gravity model. In its nonper- 
turbative form, the spectral action is defined in terms of the spectrum of the 
Dirac operator, on a spin manifold or more generally on a noncommutative 
space (a spectral triple), as the functional Tr(/(L'/A)), where / is a smooth 
test function and A is an energy scale that makes D/A dimensionless. 

The reason why this can be regarded as an action functional for gravity 
(or gravity coupled to matter in the noncommutative case) lies in the fact 
that, for large energies A it has an asymptotic expansion (see [1]) of the 
form 

(2.1) TV(/(I)/A))~ /fcA'/l^r' + /(0)CD(0) + o(l), 

A:eDimSp+ 

with fk = /q°° f{v)v^~^dv and with the integrations 



D\ 



given by residues of zeta function Cd[s) = Tr(|D|~*) at the positive points 
of the dimension spectrum of the spectral triple, that is, the set of poles of 
the zeta functions. In the case of a 4-dimensional spin manifold, these, in 
turn, are expressed in terms of integrals of curvature terms. These include 
the usual Einstein-Hilbert action 



2kq 



R^d'^x 



and a cosmological term 



70 / Vgd X, 



but it also contains some additional terms, like a non-dynamical topological 
term 



TO J R*R*^d\ 



where R*R* denotes the form that represents the Pontrjagin class and inte- 
grates to a multiple of the Euler characteristic of the manifold, as well as a 
conformal gravity term 

ao j C^upa Cf^'P'^^d^x, 

which is given in terms of the Weyl curvature tensor. We do not give any 
more details here and we refer the reader to Chapter 1 of [?] for a more 
complete treatment. 

The presence of conformal gravity terms along with the Einstein-Hilbert 
and cosmological terms give then a modified gravity action functional. When 
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one considers the nonperturbative form of the spectral action, rather than 
its asymptotic expansion at large energies, one can find additional nonper- 
turbative correction terms. One of these was identified in [^, in the case of 
the 3-sphere, as a potential for a scalar field, which was interpreted in |16j 
as a potential for a cosmological slow-roll infiation scenario, and computed 
for other, non-simply connected cosmic topologies. 

3. Geometry, topology and inflation: spherical forms 

The nonperturbative spectral action for the spherical space forms Y = 
S^/T was computed recently by one of the authors [25]. It turns out that, 
although the Dirac spectra can be significantly different for different spin 
structures, the spectral action itself is independent of the choice of the spin 
structure, and it is always equal to a constant multiple of the spectral action 
for the 3-sphere S^, where the multiple is just dividing by the order of the 
group r. This is exactly what one expects by looking at the asymptotic 
expansion of the spectral action for large energies A, and the only significant 
nonperturbative effect arises in the form of a slow-roll potential, as in [5j, 

m. 



Theorem 3.1. (Teh, [25]) For all the spherical space forms Y = S^/T with 
the round metric induced from S^, and for all choices of spin structure, the 
nonperturbative spectral action on Y is equal to 

(3.1) Tr(/(Z)y/A)) = ^ (^A'PHO) - ^A/(0)) = ^TT{f{Dss /A)), 

up to order 0(A~°°). 

Correspondingly, as explained in §5 of [16], one obtains a slow-roll po- 
tential by considering the variation i— t- + cp'^ of the spectral action 
as in [5j . More precisely, one considers a Euclidean compactification of the 
4-dimensional spacetime y x M to a compact Riemannian manifold Y x 
with the compactification of size /3. One then computes the spectral 
action Tr(/i(Dy^^i/A^)) on this compactification and its variation 

(3.2) Tr(/i((D^,5, + 02)/A2)) - TT{h{Dl^s^/A')) = Vy{^), 

up to terms of order 0(A^°°), where the potential V^((/>) is given by the 
following. 

Proposition 3.2. Let Y be a spherical space form Y = S^/T with the 
induced round metric. Let a > be the radius of the sphere and /3 the 
size of the circle in the Euclidean compactification Y x . Then the 



slow-roll potential V {(p) in (3.2) is of the form 

(3.3) VricP) = 7rAWVy(|J) + ^A^/SaWy (|^), 
where 

(3.4) Vy{x) = XyVssix) and Wy(x) = Ay W53(x), 
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with 



(3.5) 



Ay 



1 



for Y = S'/r, 



and with 




(3.6) Vss{x) = / u {h{u + x) — h{u)) du and Ws3{x) = j h{u) du 
Jo Jo 



Proof. The statement follows directly from the result of Theorem 7 of [5] 



In particular, for the different spherical forms, the potential has the same 
form as that of the 3-sphere case, but it is scaled by the factor Ay = 



Notice, moreover, that in the potential Vs3{<j)) one has an overall factor 
of (Aa)'^(A/3) that multiplies the V53 term and a factor of (Aa)(A/3) that 
multiplies the W53 term. As we observed already in [16j, when one Wick 
rotates back to the Minkowskian model with the Friedmann metric, both the 
scale factor a{t) and the energy scale A{t) evolve with the expansion of the 
universe, but in such a way that A(t) ~ l/«(i) so that the product Aa ~ 1. 
In [16] we did not need to analyze the behavior of the A/3 factor, since 
we only looked at the slow-roll parameters where that factor cancels 

out. In the spectral action model of cosmology, the choice of the scale /? 
of the Euclidean compactification is an artifact of the model, which allows 
one to compute the spectral action in terms of the spectrum of the Dirac 
operator on the compact Riemannian 4-manifold Y x S^. Eventually, the 
physically significant quantities derived from the spectral action functional 
are Wick rotated back to the Minkowskian signature case. Since in its 
nonperturbative form the spectral action functional is supposed to give a 
modified gravity action functional that works at all scales, not just in the 
asymptotic expansion for large A, it seems therefore natural to set the choice 
of the length f3 in the model so that the product A/3 ~ 1 remains constant. 

Another reason for making the assumption that A/3 ~ 1 is the interpre- 
tation given in [5] of the parameter /3 in the Euclidean compactification as 
an inverse temperature. Then, up to a universal constant, that behaves like 
the inverse of an energy scale and, when rotating back to the Minkowskian 
signature, one knows that, in the expansion of the universe the scale factor 
a{t) is inversely proportional to the temperature, so that the assumption 
A/3 ~ 1 is justified. 

With this setting, the slow-roll potential one obtains in the case of the 
3-sphere is of the form 



and §5 of [T6]. 



□ 



(3.7) 



Vy{cI)) = XyVs4(I^) 



(3.8) 
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Then one has the following result for the power spectra for the various 
cosmic topology candidates given by spherical space forms. 

Proposition 3.3. LetVs,Y{k) andVty{k) denote the power spectra for the 



density fluctuations and the gravitational waves, computed as in (1.9), for 
the slow-roll potential ((/>). Then they satisfy the power law 



(3.9) 



Vt,Y{k)^ \YVt{kQ] 



^ \ ^t,s3 + -^ log{k/ko) 



where Ay = l/#r for Y = S'^/T and the spectral parameters rig ga , n^ gs , 



c^sS'^j ^''"^ computed as in (1.12) from the slow-roll parameters (1.11) 



which satisfy ey = £53, r/y = r]s3, ^y = Cs'-^- 

To see explicitly the effect on the slow-roll potential of the scaling by 
Ay, we consider the same test functions hn{x) used in |5j to approximate 
smoothly a cutoff function. These are given by 

k=0 

Figure [1] shows the graph of hn{x) when n = 20. We use this test function 
to compute the slow-roll potential using the function V(x) -|- ^W{x), after 
setting the factors Aa = 1 and A/3 = 1, and up to an overall multiplicative 
factor of vr. We then see in Figure [2] the different curves of the slow-roll 
potential for the three cases where Y = S'^/F with F the binary tetrahedral, 
binary octahedral, or binary icosahedral group, respectively given by the 
top, middle, and bottom curve. 



4. The spectral action for Bieberbach manifolds 

We now consider the case of candidate cosmic topologies that are flat 
3- manifolds. The simplest case is the flat torus T'^, which we have already 
discussed in [16]. There are then the Bieberbach manifolds, which are ob- 
tained as quotients of the torus by a finite group action. In this section we 
give an explicit computation of the nonperturbative spectral action for the 
Bieberbach manifolds (with the exception of G5 which requires a different 
technique and will be analyzed elsewhere), and in the next section we then 
derive the analog of Proposition |3.3| for the case of these flat geometries. 

Calculations of the spectral action for Bieberbach manifolds were simul- 
taneously independently obtained in [20]. 

The Dirac spectrum of Bieberbach manifolds is computed in [2T] for each 
of the six affine equivalence classes of three-dimensional orientable Bieber- 
bach manifolds, and for each possible choice of spin structure and choice of 
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Figure 1. The test function h{x) = hn{x) with n = 20. 




Figure 2. The slow-roll potentials for the binary tetrahe- 
dral, binary octahedral, and binary icosahedral cases. 



flat metric. These classes are labeled Gl through G6, with Gl simply being 
the flat 3-torus. 

In general, the Dirac spectrum for each space depends on the choice of 
spin structure. However, as in the case of the spherical manifolds, we show 
here that the nonperturbative spectral action is independent of the spin 
structure. 

We follow the notation of [2T] , according to which the different possibilities 
for the Dirac spectra are indicated by a letter (e.g. G2{a)). Note that it is 
possible for several spin structures to yield the same Dirac spectrum. 
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The nonperturbative spectral action for Gl was computed in |il6j . We 
recall here the result for that case and then we restrict our discussion to the 
spaces G2 through GQ. 

4.1. The structure of Dirac spectra of Bieberbach manifold. The 

spectrum of the Bieberbach manifolds generally consists of a symmetric com- 
ponent and an asymmetric component as computed in [21]. The symmetric 
components are parametrized by subsets I d 1?, such that the eigenvalues 
are given by some formula Xx, x G I, and the multiplicity of each eigenvalue, 
A, is some constant times the number of x £ I such that A = A^;. In the case 
of G2, G4, G5, G6 the constant is 1, while in the G3 case the constant is 2. 

The approach we use here to compute the spectral action nonperturba- 
tively consists of using the symmetries of A^ as a function of x G / to almost 
cover all of the points in and then apply the Poisson summation formula 
as used in [5]. By "almost cover", it is meant that it is perfectly accept- 
able if two-, one-, or zero-dimensional lattices through the origin are covered 
multiple times, or not at all. 

The asymmetric component of the spectrum appears only some of the 
time. The appearance of the asymmetric component depends on the choice 
of spin structure. For those cases where it appears, the eigenvalues in the 
asymmetric component consist of the set 

B= S^27r^{kfi + c)\fiezY 

where c is a constant depending on the spin structure, and k is given in the 
following table: 



Bieberbach manifold 


k 


G2 


2 


G3 


3 


G4 


4 


G5 


6 



For no choice of spin structure does G6 have an asymmetric component 
to its spectrum. Each of the eigenvalues in B has multiplicity 2. Using the 
Poisson summation formula as in [5], we see that the asymmetric component 
of the spectrum contributes to the spectral action 

(4.1) ^ / f{u')du. 

The approach described here is effective for computing the nonperturba- 
tive spectral action for the manifolds labeled in [2T] as G2, G3, G4, G6, but 
not for G5. Therefore, we do not consider the G5 case in this paper: it will 
be discussed elsewhere. 
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4.2. Recalling the torus case. We gave in Theorem 8.1 of [16] the exphcit 
computation of the non-perturbative spectral action for the torus. We recall 
here the statement for later use. 

Theorem 4.1. Let = M^/Z^ be the flat torus with an arbitrary choice of 
spin structure. The nonperturbative spectral action is of the form 

(4.2) Tr(/(L)VA^)) = ^ / f{u^ + + w^)dudvdw, 
up to terms of order 0{A^°°). 

4.3. The spectral action for G2. The Bieberbach manifold G2 is the one 
that is described as "half-turn space" in the cosmic topology setting in [22j, 
because the identifications of the faces of the fundamental domain is achieved 
by introducing a vr-rotation about the z-axis. It is obtained by considering 
a lattice with basis ai = {0,0, H), 02 = (L,0,0), and 03 = {T,S,0), with 
H,L,S e M+ and T e M, and then taking the quotient Y = MVG'2 of by 
the group G2 generated by the commuting translations ti along these basis 
vectors a, and an additional generator a with relations 

(4.3) = ti, at20t~^ = f^"^, atsa"^ = t^^ . 

Like the torus T^, the Bieberbach manifold G2 has eight different spin 
structures, parameterized by three signs 5i = ±1, see Theorem 3.3 of [21j . 
Correspondingly, as shown in Theorem 5.7 of [21j, there are four different 
Dirac spectra, denoted (a), (6), (c), and {d), respectively associated to the 
the spin structures 





61 


S2 


'^3 


(a) 


±1 


1 


1 


ih) 


±1 


-1 


1 


(c) 


±1 


1 


-1 


(d) 


±1 


-1 


-1 



We give the computation of the nonperturbative spectral action separately 
for each different spectrum and we will see that the result is independent of 
the spin structure and always a multiple of the spectral action of the torus. 

4.3.1. The case of G2{a). In this first case, we go through the computation 
in full detail. The symmetric component of the spectrum is given by the 
data (L21J) 

/ = {{k,l,m)\k,l,m G Z,m > 1} U {{k,l,m)\k,l G Z,/ > l,m = 0} 
= ±2vry^(^ + ^)2 + i,Z2 + i_(^_p)2, 

We make the assumption that T = L. Set p = m — I. Then we have 
equivalently: 

/ = {(A;, l,p)\k, l,p£Z,p> -1}U {{k, l,p)\k, I £Z,l>l,p = -1} =: hU h 
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Theorem 4.2. Let G2{a) be the Bieberbach manifold R^/G2, with T = L 
and with a spin structure with 5i = {±1, 1, 1}. The nonperturbative spectral 
action of the manifold G2{a) is of the form 

(4.4) Tr(/(DVA^)) = HSL ( —\ [ f{u^ + v'^ + w^)dudvdw, 
up to terms of order 0{A~°^). 



Proof. We compute the contribution to the spectral action due to Ii . Since 

^klp 



Xtj is invariant under the transformation / i— t- — / and p i— )• — p, we see that 



^/(A^A') = 2^f{Xl,^/A') + fiXl,/A'). 

Z3 h P=-l 

The decomposition of Z'^ used to compute this contribution to the spectral 
action is displayed in figure |4| Applying the Poisson summation formula we 
get a contribution to the spectral action of 



plus possible terms of order 0(A~°°). 

As for I2 we again use the fact that the spectrum is invariant under the 
transformation I 1— )• — /, pi—)- —p to see that 

J]/(A2,,(_,)/A2) = 2^/(Ai,,/A2)+ ^ f{\ljA'). 

h p=l=0 

The decomposition for this contribution to the spectral action is displayed 
in figure [Sj We get a contribution to the spectral action of 



fin' 



plus possible terms of order 0(A 



When we include the contribution (4.1 ) due to the asymmetric component 
we see that the spectral action of the space G2-(a) is equal to 

TrfiD^/A') = HSL ( [ f{u' + v'^ + w')dudvdw 

again up to possible terms of order 0(A^°°). □ 
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Figure 4. Lattice decomposition for the I2 contribution to 
the spectral action of G2{a) Two regions and the set I = 0. 
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4.3.2. The case of G2{b) and G2{d). The spectra of G2{b) and G2{d) have 
no asymmetric component. The symmetric component is given by 

/ = {(A;,/,m)|A;,/,m G Z, / > 0} 
Let us once again assume that T = L. 

Theorem 4.3. Let G2{b) and G2{d) be the Bieberbach manifolds 1^3/^2, 
with T = L and with a spin structure with 5i = {±1,-1,1} and Si = 
{±1,-1,-1}, respectively. The nonperturbative spectral action of the man- 
ifolds G2{b) and G2[d) is again of the form 

(4.5) Tr(/(DVA^)) = HSL ( —\ [ f{u^ + + w'^)dudvdw, 
up to terms of order 0(A^°°). 

Proof. With the assumption that T = L and letting p = m — I, we can 
describe the spectrum equivalently by 

I = {{k,l,p)\k,l,p£Z,l>0} 



Using the symmetry 

l^ -l-l, 

we cover exactly, (see figure [s]) and we obtain the spectral action 
Tr{f{D^/A^)) = HSl( [ f{u^ + v^ + w^)dudvdw + 0{A-°^). 



□ 



4.3.3. The case of G2(c). In this case, the symmetric component of the 
spectrum is given by 

I = {(k,l, m)\k, Z, m G Z, m > 0} 



= ±2vrV ^(A: ± + + + 1/2) - ^-ly. 



1 l.o 1 ,o 1 

Again, we assume T = L. 

Theorem 4.4. Let G2{c) be the Bieberbach manifolds M.^/G2, with T = L 
and with a spin structure with 6i = {±1, 1, — !}• The nonperturbative spectral 
action of the manifold G2[c) is again of the form 

(4.6) Tr(/(L»VA^)) = HSL ( —\ f f{u^ + v"^ + w^)dudvdw, 
up to terms of order 0(A^°°). 
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1 
1 









Figure 5. Lattice decomposition for G2{b), {d) computa- 
tion. Two regions. 

Proof. If we substitute p = m — I, we see that we may equivalently express 
the symmetric component with 

I = {{kJ,p)\kJ,p£Z,p>-l} 

Using the symmetry 

I —I p 1 — p, 

we cover exactly (see figure [6|, and so the spectral action is again given 
by 

Trf{D^/A^) = HSl( f f{u^ + + w^)dudvdw + 0{A-^). 

□ 

4.4. The spectral action for G3. The Bieberbach manifold G3 is the 
one that, in the cosmic topology setting of [22] is described as the "third- 
turn space". One considers the hexagonal lattice generated by vectors 
ai = {0,0, H), 02 = {L,0,0) and as = {-^L,^L,0), for H and L in 
Ml, and one then takes the quotient of by the group G3 generated by 
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P 




1 



Figure 6. Lattice decomposition for G2(c) computation. 
Two regions. 

commuting translations along the vectors and an additional generator 
a with relations 

(4.7) = ti, at2a^^ = t^, at^a^^ = t2^t^^ . 

This has the effect of producing an identification of the faces of the funda- 
mental domain with a turn by an angle of 27r/3 about the z-axis, hence the 
"third-turn space" terminology. 

As shown in Theorem 3.3 of [21j, the Bieberbach manifold G3 has two 
different spin structures, parameterized by one sign 6i = ±1. It is then 
shown in Theorem 5.7 of [21] that these two spin structures have different 
Dirac spectra, which are denoted as G3{a) and G3{b). We compute below 
the nonperturbative spectral action in both cases and we show that, despite 
the spectra being different, they give the same result for the nonperturbative 
spectral action, which is again a multiple of the action for the torus. 

4.4.1. The case of G2>{a) and G3(6). The symmetric component of the spec- 
trum is given by 

(4.8) I = {{k,l,m)\k,l,meZ,l> l,m = 0, . . . ,1 - I}, 

(4-9) A±„ = ±27r^ ^{k + cY + -1/2 + _!_(/ _ 2m)\ 

with c = 1/2 for the spin structure (a) and c = for the spin structure ih). 
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The manifold G3 is unusual in that the multipHcity of A^^^ is equal to 
twice the number of elements in I which map to it. 

Theorem 4.5. On the manifold G3 with an arbitrary choice of spin struc- 
ture, the non-perturbative spectral action is given by 

(4.10) Tr(/(L>VA^)) = ^ (—^ HL^ [ f (u^ + + t^) dudvdt 
plus possible terms of order 0(A^°°). 

Proof. Notice that A^^^ is invariant under the linear transformations R, S, T, 
given by 

R{1) = -I 
R{m) = —m 

S{1) = m 
S{m) = I 

T{1) = l-m 
T{m) = —m 

Let / = {(A:, /, m)\k, l,m G Z,Z > 2,m = 1,.. . ,/ — 1}. 
Then we may decompose as (see figure [?]) 

(4.11) = lU R{I) U S{I) U RS{I) U T{f) U RT{f) U {I = m}. 
Therefore, we have 

5^/(Ai,^/A2) = 4j^/(Ai,^/A2) 



+ 2 E/(^fc^™M')- E fi^lim/^' 

\ I m=0, 1>1 

l=m 

= 6j;/(AL^/A2)- J]/(Ai,^/A2) 

/ m=0 
m=0, (=0 l=m 

Y^fixlu^^) = I (T^f^^lim/^') + E fi^liJ^')] 

I \ Z3 m=0 / 

E /(AL™/A^)-E/(^i^-/^')) 

Ym=0, /=0 l=m J 
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Therefore the symmetric component of the spectrum contributes to the 
spectral action 



+ 



A 

-) 



HL' I f{u'^ + + l{v-2w) 



HL 



fin' + 3^^) 



/ A \ 



V2vr J 



H fin 



- HL 



_ 4 
~ 6 
+ 0(A 



A 
2^ 



fin' 



2wy 



A 

2^' 



H fin') 



Combining this with the asymmetric contribution (4.1), we see that the 
spectral action of spaces G3{a) and G3{b) is equal to 

3 



A 
2^ 



f 



u' + v' + ^{v 



2wy dudvdw + 0{k' 



Now, if one makes the change of variables {u,v,w) i— t- (u,v,t), where 

2w — V 

t 



V3 



then the spectral action becomes 



1 

7! 



-V 

2tt J 



HL' f f{u' + v' + t') dudvdt + 0{Ar 

JR3 



□ 



Notice that, a priori, one might have expected a possibly different result 
in this case, because the Bieberbach manifold is obtained starting from a 
hexagonal lattice rather than the square lattice, but up to a simple change 
of variables in the integral, this gives again the same result, up to a multi- 
plicative constant, as in the case of the standard flat torus. 

4.5. The spectral action for G4. The Bieberbach manifold GA is referred 
to in [22] as the "quarter-turn space" . It is obtained by considering a lattice 
generated by the vectors oi = (0,0, -ff), a2 = (-L,0,0), and as = (0,L,0), 
with H,L > 0, and taking the quotient of by the group G4 generated 
by the commuting translations U along the vectors and an additional 
generator a with the relations 



at2a 



This produces an identification of the sides of a fundamental domain with 
a rotation by an angle of 7r/2 about the z-axis. Theorem 3.3 of [21] shows 
that the manifold G4 has four different spin structures parameterized by two 
signs di = ±1. There are correspondingly two different forms of the Dirac 
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m 




1 



Figure 7. Lattice decomposition for G3 computation. Six 
regions and the set I = m. The dashed hnes indicate one of 
the boundary hnes which define the region / along with its 
images under the symmetries of Xkim- The other boundary 
hne of / overlaps with the boundary of /. 

spectrum, as shown in Theorem 5.7 of [21], one for 5i = {±1, 1}, the other 
for 5i = {±1, -1}, denoted by GA{a) and G(4)6. 

Again the nonperturbative spectral action is independent of the spin 
structure and equal in both cases to the same multiple of the spectral action 
for the torus. 

4.5.1. The case ofG4:{a). 

Theorem 4.6. On the manifold G4 with a spin structure (a) with 5i = 
{±1, 1}, the non-perturbative spectral action is given by 



(4.13) 




plus possible terms of order 0(A °°). 

Proof. The symmetric component of the spectrum is given by 
I = {{k,l,m)\k,l,m e Z, / > 1, m = 0, . . . , 2/ - 1} 



= ±2vry^(fe + l)2 + -l(/2 + (^_/))2, 
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1 



Figure 8. Lattice decomposition for G4(a) computation. 
Four regions. 

First, we make the change of variables p = m—l. Then we use the symmetries 

I ^ p p ^ I 
I ^ p p^ —I 



to cover all of except for the one-dimensional lattice {{k,l,p)\l = p = 0}. 
This decomposition is depicted in figure [8} In the figure one sees that the 
points I = p such that / < are covered twice, and the points I = p such that 
/ > are not covered at all, but via the transformation i— )• —{l,p), this 
is the same as covering each of the points I = p, I ^ once. Observations 
like this will be suppressed in the sequel. Then we see that the contribution 
from the symmetric component of the spectrum to the spectral action is 

(4.14) HL^ I^ J{u^ + v^ + w^)dudvdw-^(^^^H lj{u')du, 

up to terms of order 0(A^°°). Combining this with the asymmetric compo- 



nent, we find that the spectral action is given by (4.13). □ 



4.5.2. The case of GA{b). In this case there is no asymmetric component in 
the spectrum. The symmetric component is given by the data 

/ = {(fc,/,m)|A;,/,m G Z, / > l,m = 0,. . . ,2/ - 2} 
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Figure 9. Lattice decomposition for G4(6) computation. 
Four regions. 



±27rW ^(fc + Ir + ^{{1 - 1/2)2 + (m-l + l/2))2. 



We again obtain the same expression as in the G4(a) case for the spectral 
action. 

Theorem 4.7. On the manifold GA with a spin structure [h) with 5i = 
{±1,-1}, the non-perturbative spectral action is also given by 

(4.15) Tv{f{D^/A^)) = - ( HL^ ! f{u^ + + w^)dudvdw 

2 V27ry Jm3 

up to possible terms of order 0{A^°^). 

Proof. We make the change of variables p = m — I. Using the symmetries 

l^l-l 
I p p I 
I ^ p p 1-^ 1 — I, 



we can exactly cover all of , as shown in figure [o] and so the spectral action 
has the expression (4.15). □ 



Remark 4.8. The technique we use here to sum over the spectrum to 
compute the non-perturbative spectral action does not appear to work in 
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the case of the Bieberbach manifold G5, which is the "sixth-turn space" 
described from the cosmic topology point of view in [22,] , namely the quotient 
of by the group G5 generated by commuting translations tj along the 
vectors ai = {0,0, H), 02 = (-^,0,0) and 03 = (^L, ^L,0), H,L> 0, and 
an additional generator a with = ti, at2a~^ = ts and at^a'^ = t2^t3, 
which produces an identification of the faces of the fundamental domain with 
a 7r/3-turn about the z-axis. This case will therefore be analyzed elsewhere, 
but it is reasonable to expect that it will also give a multiple of the spectral 
action of the torus, with a proportionality factor of ffL^/(4\/3). 

4.6. The spectral action for G6. We analyze here the last remaining case 
of compact orientable Bieberbach manifold G6, the Hantzsche-Wendt space, 
according to the terminology followed in p2j . This is the quotient of by 
the group G6 obtained as follows. One considers the lattice generated by 
vectors ai = {0,0, H), 02 = (-L, 0,0), and 03 = (0,5,0), with H,L,S > 0, 
and the group generated by commuting translations ti along these vectors, 
together with additional generators a, f3, and 7 with the relations 



(4.16) 





= tu 


at2a ^ 


= t2\ 


at^a ^ = ^3 ^ 




= t2, 




= t-,\ 






= t3, 












7/3a = 


- hh. 





This gives an identification of the faces of the fundamental domain with a 
twist by an angle of tt along each of the three coordinate axes. 

According to Theorems 3.3 and 5.7 of [21], the manifold G6 has four 
different spin structures parameterized by three signs 6i = ± subject to the 
constraint 616263 = 1, but all of them yield the same Dirac spectrum, which 
has the following form. 

The manifold G6 also has no asymmetric component to its spectrum, 
while the symmetric component is given by 

/ = {{k,l,m)\k,l,m e Z,l>0,k>0} 
We then obtain the following result. 

Theorem 4.9. The Bieberbach manifold G6 with an arbitrary choice of spin 
structure has nonperturbative spectral action of the form 

(4.17) Trf{D^/A^) = - ( HLS [ f{u^ + v'^ + w'^)dudvdw 

2 V2vry Jk3 

up to terms of order 0{A^°°). 
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Figure 10. Lattice decomposition for G6 computation. 
Four regions. 



Proof. Using the three transformations 

k -k - I, 
-l-l, 

k ^ —k — I / H> — ? — 1, 



one exactly covers as seen in figure llO^ and so we see that the nonper- 
turbative spectral action is given by (4.17). □ 



5. Geometry, topology and inflation: flat manifolds 



As shown in Theorem 8.3 of |T6], on a flat torus of sides £ = 1 the slow 
roll potential is of the form 



with Vx3 (x) given by 

roc 

(5.1) Vt3(x) = / u {h{u + x) - h{u)) du, 

Jo 

as in (3.6 ). 
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-3, 



Proposition 5.1. Let Y he a Bieberbach manifold Y = T /T with the 
induced flat metric. Then the slow-roll potential Vy{4>) in (3.2) is of the 
form 



(5.2) 

where 
(5.3) 



Vy(x) = Ay Vt3(x) 



with Vx'iix) as in (5.1) and the factor Ay given by 

T = G2 



(5.4) 



Ay 



/ HSL 

2 

HI? 
2^3 

HI? 
4 

HLS 



r 
r 
r 



G3 
G4 
G6 



Proof. The result follows directly from the theorems proved in §4.1| above, 
which show that the nonperturbative spectral action for y is a multiple of 
the spectral action for with proportionality factor given by Ay as in (5.4). 
The potentials 



TV(M(I^L5i -Tr(/i(D^,5i/A2)) = Vy{^) 

are then related by the same proportionality factor Ay. 



□ 



We the obtain the following analog of Proposition 3.3 in the flat case. 

Proposition 5.2. LetVgyik) andVty{k) denote the power spectra for the 
density fluctuations and the gravitational waves, computed as in (1.9), for 
the slow-roll potential Vy ((/>). Then they satisfy the power law 



(5.5) 



Vs,Y{k) 



'Pt,Y{k) 



XrVsiko) 
XyVtiko) 



l-n 



■ log(fe/fco) 



;.t3" 



^ log(fc/fco) 



where Ay is as in (5.4), for Y = T /V a Bieberbach manifold and the 
spectral parameters Ug rpijU^rs, ctgrpa, a^rp3 are computed as in (1.12) from 
the slow-roll parameters (1.11), which satisfy ey = e^a, rjy = rjrps, ^y = ^j^s. 

If one assumes that each of the characteristic sizes involved, H, L, S would 
be comparable to A~^, after Wick rotating from Euclidean to Lorentzian 
signature, as in the expansion scale A(t) ~ l/a{t), one would then obtain 
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Figure 11. The slow-roll potentials for the G2 case (top 
curve) , the G3 case (middle curve) , and the G4 and GQ cases 
(lower curve). 



proportionality factors that are simply of the form 



(5.6) 



Ay 



1 

2 

1 

2^3 

1 
4 

1 
4 



r 
r 
r 
r 



G2 
G3 
G4 
G6 



Assuming then that A/3 = 1, and using the same test function hn{x) with 
n = 20 as in Figure [T] we then obtain different curves as in Figure [TT] for the 
G2 case (top curve), G3 case (middle curve), and for the G4 and G6 cases 
(bottom curve). 



6. CONCLUSIONS: Inflation potential, power spectra, and cosmic 

TOPOLOGIES 

We have seen in this paper that, in a modified gravity model based on 
the non-perturbative spectral action functional, different cosmic topologies, 
either given by spherical space forms or by flat Bieberbach manifolds, leave 
a signature that can distinguish between the different topologies in the form 
of the slow roll inflation potential that is obtained from the variation of the 
spectral action functional. The amplitude of the potential, and therefore the 
amplitude of the corresponding power spectra for density perturbations and 
gravitational waves (scalar and tensor perturbations), differs by a factor that 
depends on the topology, while the slow-roll parameters only detect a differ- 
ence between the spherical and flat cases. As one knows from [13j , [23j, [24J, 
both the slow-roll parameters and the amplitude of the power spectra are 
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constrained by cosmological information, so in this kind of modified gravity 
model, one in principle obtains a way to constrain the topology of the uni- 
verse based on the slow-roll inflation potential, on the slow-roll parameters 
and on the power spectra for density perturbations and gravitational waves. 
The factors Xy that correct the amplitudes depending on the topology are 
given by the following table. 



Y spherical 


Ay 


Y flat 


Ay 


sphere 


1 


flat torus 


1 


lens 


1 

N 


G2{a){b){c){d) 


HSL 
2 


binary dihedral iN 


1 
AN 


G3{a){b) 


2^/3 


binary tetrahedral 


1 

24 


G4(a)(6) 


4 


binary octahedral 


1 

48 


G5 


? 


binary icosahedral 


1 
120 


G6 


HLS 
4 



Notice that some ambiguities remain: the form of the potential and the 
value of the scale factor A alone do not distinguish, for instance, between 
a lens space with N = 24, a binary dihedral quotient with = 6 and the 
binary tetrahedral quotient, or between the Poincare dodecahedral space 
(the binary icosahedral quotient), a lens space of order N = 120 and a 
binary dihedral quotient with = 30. At this point we do not know 
whether more refined information can be extracted from the spectral ac- 
tion that can further distinguish between these cases, but we expect that, 
when taking into account a more sophisticated version of the spectral ac- 
tion model, where gravity is coupled to matter by the presence of additional 
(non-commutative) small extra-dimensions (as in [H], |7]), one may be able 
to distinguish further. In fact, instead of a trivial product X x F, one can 
include the non-commutative space F using a topologically non-trivial fibra- 
tion over the 4-dimensional spacetime X and this allows for a more refines 
range of proportionality factors Ay. We will discuss this in another paper. 
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